Experimental realization of plaquette resonating valence bond states 
with ultracold atoms in optical superlattices 
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The concept of valence bond resonance plays a fundamental role in the theory of the chemical bond 
and is believed to lie at the heart of many-body quantum physical phenomena. Here we show direct 
experimental evidence of a time-resolved valence bond quantum resonance with ultracold bosonic 
atoms in an optical lattice. By means of a superlattice structure we create a three-dimensional 
array of independent four-site plaquettes, which we can fully control and manipulate in parallel. 
Moreover, we show how small-scale plaquette resonating valence bond (RVB) states with s- and 
(i-wave symmetry can be created and characterized. We anticipate our findings to open the path 
towards the creation and analysis of many-body RVB states in ultracold atomic gases. 

PACS numbers: 03.75.Lm, 03.65.Xp, 75.10. Jm, 75.10.Kt 



In his theory of the chemical bond, Pauling developed 
the concept of quantum resonance: a quantum superpo- 
sition of resonant structures with different arrangements 
of the valence bonds [1]. Such resonant states are essen- 
tial to explain the chemical properties of certain organic 
molecules like benzene [2] . In the context of high temper- 
ature superconductivity, Anderson extended Pauling's 
notion to a macroscopic level, by proposing that elec- 
trons in Mott insulating solid state materials could form 
resonating valence bond (RVB) states [3, 4]. In a Mott in- 
sulating phase, electrons are localized to individual atoms 
or molecules, and the fluctuations in the charge (density) 
degree of freedom are strongly suppressed. The physics 
is dictated by the remaining spins, which interact via su- 
perexchange interactions. Under certain conditions, the 
localized spins are expected to evade local order and con- 
tinue to fluctuate down to zero temperature, forming a 
coherent superposition of many different arrangements 
in which the spins are paired up into singlets or valence 
bonds. 

Ultracold atomic gases in optical lattices [5-7] and 
other quantum optical systems are promising candidates 
for the quantum simulation of RVB states [8] . Their real- 
ization would allow one to gain valuable insight into the 
entanglement properties of these states as well as to an- 
swer fundamental questions in condensed matter physics 
like their stability under specific Hamiltonians such as 
the Hubbard model, or to test their exotic superconduct- 
ing properties upon doping [9, 10]. In this work we cre- 
ate an array of small-scale versions of Pauling-like RVB 
states in four-site plaquettes and study their basic physi- 
cal properties. Our techniques can be directly generalized 
to a gas of fermionic atoms, for which one expects that 
the adiabatic connection of such isolated plaquette RVB 
states could lead to the creation of a macroscopic d-wave 




FIG. 1. Schematics of a single plaquette and energy levels 
at half filling, (a) Scheme of the lattice potential in the x^y 
plane, created by a pair of bichromatic optical lattices. The 
elementary cell is made of four wells arranged in a square 
configuration, (b) Energy levels of four atoms on a plaquette 
in a Mott insulating state at half filling, with superexchange 
spin couplings along x (y) denoted by Jx{Jy)- For any ratio 
Jx/Jy, the highest energy state is a total spin-| singlet. In 
the case of Jx/Jy = 0, it corresponds to the valence bond 
state 100)5 whereas for Jx/Jy = 1 it is the s-wave RVB state 
|$+). The other total singlet for Jx — Jy, lower in energy, 
is the d-wave RVB state |$-) = (c) Symbols used for 

a singlet bond and for the s-wave and cZ-wave plaquette RVB 
states. 

superfiuid state [9-11]. 

Let us consider an ultracold gas of bosonic atoms in 
two internal states, loaded into a two-dimensional su- 
perlattice structure whose elementary cell is a plaque- 
tte made out of four wells arranged in a square pattern 



[Fig. 1(a)]. In the regime in which the tunnehng amph- 
tude between adjacent plaquettes is strongly suppressed, 
the system can be regarded as a cohection of independent 
rephcas of a single plaquette, the object of our study. At 
half filling, and when the on-site interaction U dominates 
over the tunneling amplitude t between wells in a plaque- 
tte, atoms are site localized, one per site, and the physics 
is governed by the remaining four effective | -spins, which 
interact with their next neighbors via a ferromagnetic 
Heisenberg interaction JSi-Sj^ with J = —4£^/U [12-16]. 

To gain insight into the RVB states on a plaquette, it is 
convenient to write the Heisenberg interaction in terms of 
the swap operator Xij = 2Si-Sj + 1/2, a unitary operator 
that exchanges the states of the spins on the sites i and 
j. The plaquette Hamiltonian then takes the form [Al]: 



JyXy^ 



(1) 



where X^(^y^^ involves exchanges of two spins along an x 

(^)-bond: X, = {Xab^Xcd)/2, Xy = {Xad^Xbc)/^, 
with A^B^C^D labeling the four sites of the plaquette 
[Fig. 1(a)]. From now on, we consider solely the subspace 
of total spin zero, where all spins are part of a singlet 
state or valence bond. This subspace is generated by 
two states, which correspond to arrangements in either 
vertical |00) or horizontal bonds |g) [Fig. 1(b)]. 

Within this subspace and for identical superexchange 
couplings Jx = Jy = J, the Hamiltonian of Eq. (1) re- 
duces to = —JXxy^ where Xxy = {^ac + Xbd)/'^ 
swaps two spins along a diagonal. As can directly be 
seen, this diagonal exchange is equivalent to a 90-degree 
rotation of the plaquette and converts the state |00) into 
|§) and vice- versa, giving rise to a resonance. The eigen- 
states are then coherent superpositions of the form: 



± 



These minimum instances of RVB states exhibit no local 
magnetic order, and can not be distinguished from each 
other by measuring single-site spin observables. How- 
ever they are distinct with respect to an exchange of two 
spins along a diagonal: the 5- wave RVB state |$+) is 
symmetric; the (i-wave RVB state \^-) is antisymmetric, 
owing to its singlet structure along the diagonals of the 
plaquette, | <!>_) = 

Our experiments began with a quasi-pure Bose- 
Einstein condensate of about 5 x 10^ ^^Rb atoms in the 
Zeeman state |F = 1, ttif = — 1). The atoms were loaded 
into a tetragonal optical lattice potential, formed by three 
mutually orthogonal standing waves with wavelengths 
As = 767 nm ("short lattices") along x and and = 
844 nm along z. Two additional standing waves with 
wavelengths of A^ = 1534 nm ("long lattices") that were 
superimposed with the short lattices [18] along x and y 
were then used to realize a three-dimensional periodic 
potential whose elementary cell is a plaquette [Fig. 1(a)]. 
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FIG. 2. Initial state preparation and valence bond oscilla- 
tions, (a) Schematics of the preparation of an array of va- 
lence bond states |00) from a unit-filling Mott insulator. (b) 
Schematics of the valence bond oscillation: starting from 
100)? we switch on identical superexchange couplings along 
X and ^, leading to a coherent oscillation between |00) and 
(c,d) Fraction of band excitations rigxc (c) and STO 
amplitude ^sto (^) ^ function of the hold time r, with 
Jx — Jy = h X 120(10) Hz. (e) Frequency of the valence bond 
oscillation as a function of Jx/ Jy (f) Ratio Ay^/ (Avb+^vb) 
as a function of Jx/Jy, where Ay^ is the initial amplitude of 
the valence bond oscillations as shown in (c). The solid lines 
in (e) and (f) are calculated from Eq. (1). The horizontal er- 
ror bars represent the uncertainties in lattice depths and the 
vertical ones represent the la uncertainties of the fits to the 
STO traces. 



The final lattice depths were chosen to access the Mott in- 
sulating regime with at most one atom per lattice site for 
our total particle number. We then employed a sequence 
of site merging, spin changing collision (SCC) [19] and 
singlet-triplet oscillation (STO) operations [21, Al, A2] 
on plaquettes [see Fig. 2(a)] in order to create the initial 
state 1 00) out of the atomic spin states \F = l,mi? = —1) 
and \F = I^tuf = +1) [22]. In total we operate in paral- 
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lei over about 10^ identical plaquettes with unit atom fill- 
ing. Lattice depths of V^i = Vyi = 35 and = 40 
ensure negligible atom tunneling between plaquettes [23] . 

To directly observe the valence bond resonance, the 
initial state |00) was evolved under the Hamiltonian of 
Eq. (1) with identical superexchange couplings along x 
and y. To this aim we ramped down the short-lattice 
depths in 200 fis to Vxs = Vys = 12 Ef., resulting in equal 
couplings Jx — Jy = h X 120(10) Hz and a suppression of 
first order tunneling as t/U 1/8. Since X^y = 1, the 
evolved quantum state at time r is 

l*(T))=cos^|00)-isin^|§), (2) 
oscillating between the states |00) and |g) with frequency 

^ = 2Jx{y)/h. 

To characterize this state evolution, we measured the 
projections onto the two valence bond states: Cx = 
|(^(r)|g)|^ and Cy = |(^(r)|00)|^ which are expected 
to show oscillations of amplitude 3/4, since (g|00) = 1/2. 
Within the subspace of total singlets, the observable Cx 
can be obtained either by measuring the fraction of band 
excitations n^^^ = Cx/2 after merging pairs of wells along 
the X direction, or by measuring the amplitude Agr^Q of 
STO [21, Al, A2] induced by a magnetic- field gradient 
along X [22]. As shown in Fig. 2(c) and (d), we indeed 
observed a coherent evolution of both n^^^ and ^gxo- 
This dynamics corresponds to anti-correlated oscillations 
of the projections Cx and Cy that reveal the periodic 
swapping of the valence-bond direction. The measured 
oscillation frequency a;/27r = 250(10) Hz is compatible 
with twice the value of the superexchange couplings, in 
agreement with Eq. (2). While the damping of the va- 
lence bond oscillation (1/e decay time of 6(1) ms) could 
be attributed to inhomogeneities of the different plaque- 
tte parameters across the atomic sample, the slow overall 
increase of n^^^ and n^^^ could be caused by decoherence 
within a plaquette. We provide further evidence of the 
valence-bond dynamics governed by superexchange in- 
teractions by studying the dynamics for anisotropic cou- 
plings Jx 7^ Jy As shown in Fig. 2(e) and (f), the mea- 
sured oscillation frequencies and amplitudes as a function 
of Jx/Jy agree well with the values predicted from the 
Hamiltonian dynamics of Eq. (1). Site-resolved popula- 
tion measurements were used to check that throughout 
the evolution the four plaquette sites remained equally 
populated [18, 22]. In the absence of residual magnetic 
field gradients we expect the atoms to remain in the sin- 
glet subspace 5* = 0. This was checked by holding singlet 
atom pairs after the initial state preparation, and observ- 
ing no conversion to triplet pairs. 

In order to create the 5- wave RVB state we made 

use of the fact that it is adiabatically connected to the 
initial state |00) [Fig. 1(b)]. To follow this adiabatic path 
we started from a situation in which Vxs = 22 and 
Vys = 12^^. For these parameters, Jx/Jy is negligible 



and 100) is an eigenstate of the Hamiltonian in Eq. (1). 
We then decreased Vxs to 12 Ef. within 5 ms using an 
exponential ramp, converting the initial state into the 5- 
wave RVB state. In order to check the adiabaticity of 
the lattice-depth ramps, we then increased the short lat- 
tice along X {y) to 22 E^ in 5 ms , transforming the RVB 
state back into a valence-bond state |00) (or respec- 
tively). By using STO we measured the singlet correla- 
tions along both directions x and y for the initial, inter- 
mediate and final states of the ramp |00) ^ |*^+) ^ |00) 
or |g) [Fig. 3(a)]. As expected, for the initial state we 
observe oscillations close to maximum amplitude only 
along y and none along x. In the intermediate state, 
the oscillation amplitudes are approximately equal, as 
expected for a non-degenerate eigenstate of the Hamil- 
tonian in Eq. (1) with symmetric couplings. After the 
second ramp, depending on whether the superexchange 
coupling was decreased along x or y^ we observe singlet 
correlations mostly along the direction of strong coupling. 
The measured amplitude of STO in the final state was 
found to be smaller than in the initial state, due to deco- 
herence in our atomic sample which occurred on a time 
scale of 30 ms in our setup. As can be seen in Fig.S2 
of the supplementary material, for a total ramp time of 
10ms (gray bar) the value of Agr^Q is 0.22(2), which is 
comparable to the STO amplitude of 0.24(2) obtained for 
the initial state |00) [see Fig. 2(d)]. 

In the RVB state the projections on the valence 

bond states are given by C^^ = C^y = 3/4. They can be 
obtained from the STO amplitudes according to Cx,y = 
1/4 + 3/2 Ag^^o [22]. By averaging the measured STO 
amplitudes around Jx = Jy^ we obtain Cx = Cy = 0.76{7) 
[Fig. 3(b)], in good agreement with the theoretical pre- 
diction. We also measured Cx^ Cy as a function of the cou- 
pling anisotropy Jx/Jy, by following the adiabatic path 
100) |^+) \^{Jx/ Jy)) with a fixed total ramp time 
of 10 ms. As shown in Fig. 3(b), the measurement results 
are in good agreement with the theoretical values in the 
adiabatic limit (solid lines) and with a model taking into 
account the finite ramp time (shaded lines). 

The d-wave RVB state |<l>_) is obtained from the state 
100) by exchanging two spins along a bond in the x direc- 
tion: 

XM) = Iss). (3) 

This unitary operation was implemented by a quantum 
evolution of the state |00) under the Hamiltonian Eq. (1) 
for = 0, yielding: 

|*(t))=cos^|00>-isin^|g2>, (4) 

with uo = 2Jx/h. For a hold time r = ir/u the initial state 
evolves into |^), characterized by C^^ = ^y = 1/4, and re- 
duced STO amplitudes ^gxo ~ ^sto = 1/8- As shown 
in Fig. 4, in that state the amplitude of STO was indeed 
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FIG. 3. Preparation of the s-wave RVB state and adiabatic 
valence bond swap, (a) Scheme of the adiabatic conversion 
100} |^+) 100) or |§), together with STO for each state. 
The STO period along x is not constant due to an increasing 
magnetic gradient field during the measurements that was 
caused by a finite response time of the coils producing the 
magnetic- field gradient along x. The solid lines are fits of the 
STO taking into account the finite rise time of the magnetic 
field gradient for the x direction, (b) Projections Cx, Cy on 
the valence bond states as a function of the ratio Jx/ Jy of su- 
perexchange couplings, measured from the STO amplitudes. 
The latter were rescaled in order to give the expected value 
of 0.5 for the valence bond states |g) and 100), using the data 
points at Jx/ Jy = 0.006(2) and Jx/ Jy = 10(3). For the point 
at Jx/ Jy — 200(50) the rate of change of the couplings was 
the largest and adiabaticity was not maintained. The hori- 
zontal error bars represent the uncertainties in lattice depths 
and the vertical ones represent the la uncertainties of the fits 
to the STO traces. The solid lines are calculated from the 
eigenstates of the Hamiltonian in Eq. (1). The shaded lines 
are calculated by modeling the experimental ramps using the 
Schrodinger equation. Their widths represent the uncertain- 
ties in the lattice calibration. 



much reduced, in our case below the noise level. However, 
the large STO amplitude along observed both in the 
initial state and after one period of evolution (r = 27r/a;), 
demonstrates the coherence of the evolution and rules out 
a reduction of contrast at r = tt/cc; due to decoherence. 
Alternatively, after preparing the |<l>_) state, we inverted 
the coupling direction by increasing in 200 /is the short- 
lattice depth along x to 22 and decreasing the one 
along y to 12 As shown in Fig. 4, we then observed 
a coherent evolution to a state with a large overlap with 
|g), according to the measured STO. 

In conclusion, we have shown direct experimental ev- 
idence of a valence-bond quantum resonance in an ar- 
ray of replicas of optical plaquettes, preparing and de- 
tecting minimum versions of RVB states. The s-wave 
and d-wave plaquette states created here could be used 
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FIG. 4. Preparation of the d-wave RVB state. Schematics 
of the experimental sequence: starting from the state |00) 
we suddenly switch on the superexchange coupling along x. 
The (i-wave RVB state |$-) is obtained at the tt time of the 
subsequent periodic evolution. Measured STO are shown at 
the 0, TT and 2n times of the evolution. For the lower right 
state, we inverted the coupling direction at the tt time. 



to encode a minimum instance of a topologically pro- 
tected qubit. When stabilized by a Hamiltonian H = 
J{Xx^Xy^Xxy), corresponding to a situation in which 
superexchange interaction takes also place along the di- 
agonal bonds, these two states form a degenerate two 
level system which is immune to local decoherence aris- 
ing, for instance, from on-site fluctuations of the external 
magnetic field. Such an arrangement could also be di- 
rectly adapted to a setting of four coupled quantum dots 
to realize protected qubits in a solid state setting [24]. 
Further extensions enabled by this work include the adi- 
abatic connection of the plaquette RVB and valence bond 
solid states, or the study of their non-equilibrium dynam- 
ics upon instantaneous coupling in quantum ladders or 
extended two-dimensional systems. Moreover, the pla- 
quette tools developed here could be used as building 
blocks for more complex protocols leading to a variety 
of topologically ordered states, like Laughlin states or 
string net condensates [25, Al]. Finally, we note that all 
presented results could also be obtained using fermions 
instead of bosons, where the singlet valence bond is the 
true ground-state of a two-spin dimer. In that case, the 
adiabatic connection of isolated RVB states could lead to 
the formation of a d-wave superfluid upon doping [9-11]. 

This work was supported by the DFG (FOR635, 
FOR801), the EU (STREP, NAMEQUAM, Marie Curie 
Fellowship to S.N.), and DARPA (OLE program). M. 
Aidelsburger was additionally supported by the Deutsche 
Telekom Stiftung. 
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APPENDIX 



FILTERING SEQUENCE 

The study presented in the main text relies on the 
loading of plaquettes at half filling, i.e. with four atoms 
in total per plaquette. Despite the preparation of the 
atomic sample in a Mott insulator state at unit filling, 
defects are relatively likely in our system. In order to iso- 
late the signal from correct configurations, we perform a 
filtering sequence that consists in transferring the atoms 
in plaquettes with incorrect fillings into different hyper- 
fine states, which are not probed in the final atom imag- 
ing [see Fig. Al(a),(b)]. We first merge pairs of sites 
along X and perform spin-changing collisions (SCC) to 
convert pairs of atoms in \F = I^tuf = 0) to the Zee- 
man states \F = l^mp = 1) and \F = l^mp = —1). A 
microwave pulse then transfers the remaining atoms in 
\F = l,mF = 0) to the state \F = 2,mF = 0). We then 
use SCC to transfer back atom pairs in |F = 1, = 0). 
With another microwave pulse we transfer the remaining 
atoms in \F = l^mp = I) due to finite SCC fidelity to 
\F = 2, mi? = 2). Then we split the sites along x. The 
rest of the sequence is identical to the one described in 
Fig. 2(a) in the main text, i.e. we merge pairs of sites 
along y and perform SCC. As shown in Fig. Al, the pla- 
quettes with one atom per site end this filtering sequence 
in the desired configuration. On the contrary, a hole in 
the initial configuration leads to the transfer of one atom 
to the F = 2 hyperfine state, and finally to a configu- 
ration with one atom per site in |F = 1, m^? = 0) before 
the final SCC. Therefore no atom is transferred to the 
state \F = l^rriF = I) that we probe at the end of the 
experiment. Similar conclusions can be obtained for the 
configurations with additional holes or particles. In to- 
tal about 10% of the atoms end the filtering sequence in 
\F = l,mF = 1) [see Fig. Al(c)]. 
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FIG. Al. Schematics of the filtering sequence for plaquettes at half filling, (a), (b) Filtering sequence applied 
to a plaquette initially filled with 4 atoms (a) and 3 atoms (b). In the case of 3 atoms per plaquette no atom is eventually 
transferred in the Zeeman state |F = 1, = 1) that we probe at the end of the experiment, (c) Evolution of the atom fraction 
in \F — l,mi? = 1) (within the F —1 manifold) during the three successive SCC used for the filtering sequence. 



MEASURING SINGLET CORRELATIONS 

The quantum states prepared in this work were probed 
by measuring the projections Cx^y on the valence bond 
states along both directions. The observable Cx can be 
measured as follows: pairs of sites are merged along x by 
decreasing the short-lattice depth along x to in 10 ms. 
For a pair of atoms in a spin-triplet state, both atoms 
are transferred to the lowest Bloch band; on the con- 
trary, due to the different parity of the quantum state, 
for a spin-singlet pair one expects one atom to occupy the 
first excited band [Al, A2]. A subsequent band-mapping 
technique allows one to measure the fraction n^^^ of band 
excitations and to infer the value of Cx = 2 71^^^. Alter- 
natively, the singlet correlation along x can be probed 
through the amplitude AgrpQ of singlet-triplet oscillations 
[A2]. We applied a magnetic-field gradient along x in 
order to induce a coherent oscillation between |c^) and 
1'=') states (here |'=') denotes a spin-triplet along the hor- 
izontal direction). An explicit calculation of the STO 
amplitude in the singlet state of highest energy gives 
Cx,y — 1/4 + 3/2 ^gxo* 

In the initial state |00), the measured STO amplitude is 
about half of the expected value Ag^pQ = 1/2. This can 
be attributed to residual excitations introduced by the 
site merging, to the residual spatial overlap after time- 
of-flight between atoms from the ground and first excited 
band, as well as to the presence of residual holes in the 
plaquette that do not contribute to STO. 

For the measurement of Cx^y as a function of Jx/ Jy [see 
Fig. 3(b) in the main text], we followed the adiabatic path 
100) \^+) ~^ \^{Jx/ Jy)) with a fixed total ramp time 
of 10 ms. The STO amplitudes were rescaled in order to 
give the expected value of 0.5 for the valence bond states 
|g) and 100), using the data points at Jx/Jy = 0.006(2) 
and Jx/Jy = 10(3). The rescaling factors were 3.2 and 
2.3 for Cx and Cy respectively. For the point at Jx/Jy = 
200(50) the rate of change of the couplings was the largest 
and adiabaticity was not maintained, as indicated by the 



theoretical model in Fig. 3(b) in the main text. 



DECOHERENCE IN THE SINGLET SUBSPACE 

We investigated the decoherence of the highest energy 
total singlet state by measuring the STO amplitude of 
the adiabatic sweep |00) ^ |*^*+) |§) as a function of 
the total ramp duration [see Fig. A2]. The decrease in 
fidelity for small ramp durations is well accounted for by 
a numerical calculation of the evolution of the quantum 
state during the ramp, according to the Hamiltonian in 
Eq. (1) in the main text. For the ramp duration of 10 ms 
used for Fig. 3 in the main text, the calculation predicts a 
AgrpQ = 0.48 [see inset Fig. A2]. The measured decrease 
of fidelity for longer ramp durations illustrates a deco- 
herence mechanism in our system whose understanding 
would require further studies. 



SITE-RESOLVED DETECTION 

To detect the atom numbers on the different sites of 
the plaquette we apply two mapping sequences along the 
X and y direction during which the populations Nq are 
transferred to different Bloch bands analog to the tech- 
nique described in Ref. [A3] for isolated double wells. 
A subsequent band-mapping technique allows us to de- 
termine the population in the Bloch bands by count- 
ing the atom numbers in different Brillouin zones. The 
colors used in Fig. A3 (a) and A3(b) show the connec- 
tion between the Brillouin zones and the correspond- 
ing lattice sites. A typical image obtained after 10 ms 
of time-of- flight is shown in Fig. A3(c). The popula- 
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FIG. A2. Decoherence in the singlet subspace. STO 

amplitude after an adiabatic sweep |00) ~^ IS)? 

a function of the total ramp duration. The solid line is the 
product of the sweep fidelity expected without decoherence 
and of an adjustable exponential decay (1/e decay time of 
27(5) ms). The gray region indicates the ramp duration for 
experiment described in the main text. Inset: STO ampli- 
tude expected without decoherence, calculated by solving the 
Schrodinger equation with the Hamiltonian shown in Eq. (1) 
in the main text. 



spin-dynamics during the oscillation. 
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FIG. A3. Site- resolved detection, (a) Schematic of the 
four-site plaquette. (b) Brillouin zones of the 2D lattice, (c) 
Typical momentum distribution obtained after 10 ms of time- 
of- flight, (d) Measured population imbalance during the va- 
lence bond oscillations for the trace shown in Fig. 2(c) in main 
text. 



tion imbalance during the valence bond oscillations is 
shown in Fig. A3(d). The population in four plaque- 
tte sites remained equally populated, proving the purely 
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